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A GENERAL SYSTEM FOR CALCULATING BURNING RATES OF PARTICLES AND IRGPS 
AND COMPARISON OF CALCULATED RATES FOR CARBON, 

BORON, MAGNESIUM, AND ISOOCTANE 
By Kenneth P. Coffin and Richard S. Brokaw 


SUMMARY 

A system with general equations has been devised for computing the 
burning rates of small particles burning as diffusion flames; the equa- 
tions account for the effects of diffusion and dissociation with a high 
degree of rigor. Two types of solutions are carried out: (l) a numeri- 

cal integration of considerable complexity, and (2) a somewhat less 
complicated and less rigorous analytical solution involving stepwise 
iteration across the temperature profile. The direct results of the 
calculation are partial pressures sis a function of temperature. Simple 
additional calculations produce the burning rate and the flame 
structure. 

Both solutions were obtained for carbon burning in air; the differ- 
ences sure slight. For boron, because of the greater number of equilibria 
involved, only the analytical solution was undertaken; a special treat- 
ment for a solid reaction product, boric aside (B 2 0 3 ), was required. A 

number of ambient temperatures and ambient oxygen concentrations were 
examined in the case of boron; as an example, three graphs rapidly yield 
burning rates for a wide range of ambient conditions. 

The general equations presented here reduce to the much simpler 
equations used by previous investigators. The simplified equations were 
also applied to boron, and yielded results in general agreement with the 
more detailed analytical solutions. Earlier results from. the simplified 
equations for isooctane and magnesium are included for comparison. The 
simplified equations appear to be sufficiently accurate for many purposes. 

For a series of substances covering a wide range' of volatility, 
relative heat-release rates are in the order.: hydrocarbon > magnesium > 

carbon “ boron. 
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introduction 

In the last few years , a number of Articles have appeared which 
treat the combustion of single fuel particle's in quiescent air. The 
treatments of these diffusion flames vary in generality and in specific 
approach (refs. 1 to 6). They involve physical and mathematical approxi- 
mations which are qualitatively acceptable, but which may introduce sig- 
nificant quantitative errors. ' 

The general success of the stagnant-jfiim treatment for the combus- 
tion of single liquid drops has tended td support the validity of the 
assumptions involved. The model assumes !a Steady-state diffusion flame 
surrounds the fuel drop. Cxygen diffused inward, fuel vaporizes and' 
diffuses outward, and they meet and react in a flame front symmetrically 
surrounding the drop; sufficient heat is [conducted inward from the reac- 
tion zone to vaporize the fuel and heat the vapor to reaction tempera- 
ture, while the remaining heat is conducted outward and, in part, is used 
to preheat the inbound oxygen. The reaction products diffuse outward. 

The assumption of a diffusion flame Wehanl mn with the over-all 
process controlled by the physical processes of heat conduction and dif- 
fusion carries the implicit assumption thht the chemical processes within 
the flame front are fast. Such assumptions are in contrast to those made 
in considering ignition delay, in which chemical processes are presumed to 
control. In the case of flame speed, as well as in quenchl ng-di stance 
and dead-space phenomena, both chemical and physical processes are pre- 
sumed to be involved. In the steady-state burning of a liquid monopro- 
pellant drop, a diffusion flame is not present, and both physical and 
chemical processes contribute to the over4all rate. 

On the basis of the diffuslcn-flame (stagnant-film) model, burning 
rates in reasonable agreement with experiment~have been calculated by 
many workers, although certain difficulties appear in the case of calcu- 
lated flame structures . The main difficulty _is that computed flame diam- 
eters are significan t ly greater than those- observed experimentally; how- 
ever, appropriate corrections involving the experimental boundary condi- 
tions result in only relatively small changes" in the burning rate. 

The general objective of this investigation was a more careful solu- 
tion of the humlng-rate problem in an effort! to determine the limita- 
tions of the simpler treatments. More specifically, this investigation 
was concerned with the calculation of the relative burning rates for 
several solid combustibles representing a wide range of volatilities. 

It seemed appropriate to make efforts I to extend the calculations 
made on liquid fuels and carbon (refs. 1 and 5) to systems of solid com- 
bustibles which do hot involve conventional fuels. The combustion of 
magnesium had been treated experimentally and a tentative approach to 
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the theoretical treatment given (ref. 6). Magnesium is a combustible 
metal of relatively low volatility (b.p., 1393° K) and as such represents 
a transition to the less-volatile metals. For the present treatment, 
boron (b.p., 2550° K) was selected as a typical high-boiling-point com- 
bustible. 

In this report, expressions for diffusion and heat c o nd u ction , as 
they apply to bu rning particles, are developed into a system of gen- 
eralized equations somewhat along the lines indicated by Spalding in his 
approximate trea tment of carbon (ref . l) . These equations are applied 
to boron and carbon, and they are reduced to the simpler equations used 
in the treatment of hydrocarbons and magnesium. Hie simplified equations 
are also applied to boron. 

Hie general equations are developed in the next section. In the 
section following that development, carbon is considered as a detailed 
specific examp le of the use of the general equations. Carbon is treated 
in two ways, in terms of an exact numerical integration and in terms of 
a simpler analytical solution involving certain average values and limit- 
ing forms of the diffusion equations. 

The boron system involves a more complex series of equilibria than 
that of the carbon system; therefore, only the analytical solution was 
undertaken. These results for boron are in turn compared with those ob- 
tained for boron by the simpler equations of the form used for isooctane 
(refs. 1 and 5) and magnesium (ref. 6). The details of Hie boron case 
are in appendix B. The general equations are reduced to the simpler 
equations for isooctane and magnesium in appendix C. Hie bur n ing rates 
and particle lifetimes of carbon, magnesium, boron, and lsooctane are com- 
pared. 


GENERAL EQUATIONS 

The development in this section deals with the application of the 
equations for the conservation of mass, the conservation of energy, dif- 
fusion, heat conduction, and chemical equilibrium to the diffusion flame 
surrounding a burning particle or drop (fig. l) . In part. It is an ex- 
tension of the early work of Spalding (ref . l) . In this section the 
equations are presented in general terns; in the following section, the 
burning of carbon is treated as a specific example of these general ex- 
pressions. Equations are designated so that numbers in this section, the 
next section, and appendixes B and C apply to corresponding equations. 

Hie general physical model Is indicated in figure 1. The diffusion 
flame symmetrically surrounding the fuel particle may be divided into 
arbitrary zones for descriptive purposes and approximate calculations; 
the boundaries of the zones are Indicated by concentric circles. 
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Convection is assumed to be absent. The ! central crosshatched portion 
represents the fuel: C represents the boundary with the ambient atmos- 
phere. B (fig. 1(a)) or the region BB* (IfigT 1(b)) represents the flame 
front, of infinitesimal or finite width, 'respectively; only di ffusion 
occurs in the finite flame front. Chemical reaction rates are assumed 
to be fast relative to heat and mass transfer. In the region AB, suffi- 
cient heat is conducted inward to vaporizje the fuel and raise it to the 
flame temperature; fuel is vaporized at t|he fuel surface A and diffuses 
to the flame front. In the regions BC of B'C, outside the flame front, 
the heat of combustion is dissipated outwjard by conduction, the reaction 
products diffuse out, and oxygen diffuses! in, with an exchange of en- 
thalpy between them. All these processes! occur in a stagnant film. The 
ambient atmosphere is a source of oxygen and a sink for heat and combus- 
tion products. 


One of the major advantages of the equations presented in this sec- 
tion is that the artificial zones described above are eliminated. This 
in turn eliminates the arbitrary assignment of flame-front boundary con- 
ditions. Of course, other assumptions arfe introduced, for example, the 
existence of equilibrium throughout the system. Further, a more rigorous 
treatment of diffusion can be employed, although difficulties exist with 
the extrapolation of diffusion constants and thermal conductivities to 
high temperatures. ! - 


The equations required may be developed 'in the following manner. 
(Carbon, as a specific example of the general equations, is treated in the 
next section; boron is treated in detail in appendix B.) X,Y,Z, . . . 
are the s atom, species present in the system; fuel or oxidant or both 
may contain several species. For example; for carbon burning in air, 
carbon, oxygen, and nitrogen are s » 3 atom species. x k i Y l i z m i • ■ • 

is the I* 111 species of the v molecular ‘species present (k^, l±, m-^, . . 
« 0, 1, 2, 3, . . .). For example, in the carbon and air case, if the 


ith 


species is the carbon dioxide molecule, C^OgNo represents C0 2 with 
2, and i hqq - 0. (Symbdls are defined in appendix A.) 


kcOg * 1 f l C0 2 


For the conservation of mass, there is a continuity equation for 
each of the s atom species; in terms of ithe total mass-transfer rate 
% " Wj m °l es per second, for ‘each of the molecular species 

across any symmetrical surface around the particle, these equations are 
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where Wf is the burning rate of the fuel in moles per second. The co- 
efficient of Wf appropriate for each equation is determined by the 
composition of the fuel molecule f ** Xj^Y^Z^. 


For systems with gaseous combustion products, there is a net out- 
ward flow of fuel Wf, because the drop is a fuel source; further, the 

net flow of oxidant and diluent gas from the ambient atmosphere is zero, 
because there is no sink in the system. The presence of some atom 
species both in the atmosphere and in the fuel (e.g., oxygen, when al- 
cohol burns in air) is provided for in the equations. Should a combus- 
tion product stick on the surface or dissolve in the fuel, a sink would 
then exist; the coefficient of Wf far the atom species involved in such 
a product would be negative and determined by the stoichiometry. 

The diffusion equations for the problem involve certain approxima- 
tions which are treated in appendix D. This treatment introduces approx- 
imations which are less severe than those made in the previous treat m e n t 
of such problems. The diffusion equation for the f 01 molecular species 
is given by 


-G ± 


Df dpj 
RT dr 


Pj ^ dp^ Wf 

pjj 2-/ HE dr " “ 4jcr 2 
c j=l 


J ?* N 2 and i j* N 2 (2a) 


(eq. (D5) in appendix D). For the case of large fractions of inert gas, 
p 4 -+ 0. Then as in equation (D7) 


r D i d Pi 
" G i " RT dr 



1 5**2 




(2b) 
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Similar equations are obtained for v - l! molecular species. If the 
entire process occurs at constant pressure, the partial pressure erf* the 
inert gas is given by 


^2 " P - i&>! P i (3) 

If no inert gas is present,, the equations (2a) must be revised in 

Pi 

terms of equations (D4) including the tent tr G v with some convenient 

reacting species designated as the v"^ 1 species. If more than one Inert 
gas is present, equations (2a) present no [difficulty; however, when the 
continuity equation Wi ne rt = 0 is used, equations (2b) indicate that 
the partial pressure of such an inert gas [is constant and equal to the 
partial pressure in the boundary condition. 

In addition, there are equilibrium equations for v - s molecular 
species (i.e., number of molecular species minus number of atom species). 
The v - s equilibrium equations are written in terms of s molecular 
species where s equals s minus the number of inert species. The s 
independent molecular species are usually Selected for convenience as those 
naninert species having the highest partial pressures (however, they must 
represent all naninert atom species) . A bjalanced chemical equation for 
the I t]l dependent species in terms of eafch (g^h) the s independent 
species can be written ! 


S : 

1 + £ n gS - 9 — (*) 

gal 

where the rig's are the stoichiometric coefficients. Then 

pi - p g ng - < 5 ) 

g»l i 

where 

s 

log JtjL - log + £ nig log Kg (6) 

g*»l i 

The K's are the equilibrium constants for, the' formation of the mole- 
cule from gaseous atoms (ref. 7), 


t 
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The amount of heat transferred by conduction is given in terms of 
the flows and enthalpies of the molecular species by 

«A 0 - £ w i% - - (?) 

Hf 0 is the enthalpy of the fuel at its initial conditions and is a 
special case of Hj_, the enthalpy of the i^ 1 species at temperature T; 
is of reference 7 and includes the chemical enthalpy. 

And so for unknowns, there are v + 1 flows (W^, Wg, . . . W^, • . . 
W v , Wj) and v pressures (p]_, p 2 , . . . Pj_, . . . p^ ) . There are, 

therefore, 2v + 1 unknowns in the system of equations. There are s 
continuity equations (eqs. (l)), v - 1 diffusion equations (eqs. (2) 
a total-pressure equation (eq. (3)), v - s equilibrium equations (eqs. 
(5)), and an energy equation (eq. (7)) for a total of 2v + l equations. 
At least in principle, the system of equations may be solved for all the 
p i and at any specified temperature. 

The elimination of Wf from the energy equation (7) and the con- 
tinuity equations (l) yields a modified energy equation and s - 1 con- 
tinuity equations. The diffusion equations are then used to eliminate 
the W's from the energy and continuity equations. The symbol 
8^ a D^/xRT (appendix E) is introduced for convenience throughout, and 
dT/dr is introduced into the continuity equations. The results are an 
energy equation and continuity equations in terms of 5^, E^, and 

dpi/dT. 

If the diffusion equations are taken in the form of equation (2a), 
the equations can be integrated numerically to produce a solution which 
is essentially exact (subject to the approximation in the transformation 
of appendix D) . Seme complications are encountered; an analytical solu- 
tion of the equations is not possible. However, a numerical solution 
has been obtained for the carbon case; the significant equations and the 
difficulties are indicated in the next section. 

If the diffusion expressions are taken in the form of equation (2b), 
the equations take the following form. Hie energy equation is 

E dPi 

8 i^L dT “ 
i 


-1 


(S) 
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and s - 1 continuity equations t- - 

(9) 

In a specific case, E£ and c^ develop ;in the algebra. 


K[ and 5^ appearing in equations :(8) and (9) are slowly varying 

functions of temperature over wide ranges of temperature,. In addition, 
8jL may vary with composition. For narrcjw temperature ranges, 8^ nwri 

H£ may be replaced by their mean values (in practice here, the 8^'s 

have been taken as constants). Then equations (8) and (9) may be inte- 
grated in closed form. 

The energy equation (eq. (8)) becomes 

"♦Swv • 

Ap l " -8-jHI 1 ^ N 2 ( i0 ) 

where AT - T - T Q and « p ± ^ T - p^jp 


form 


The continuity equations (eqs. (9)) yield s - 1 equations of the 


v —1 


Ap 


]£ c^Ae^ i ^ K [ and g a 2, 3, ... s (ll) 


8 °e 8 g i,i g 

Further, there are v - s equi lib rium equations 


pi - ^iTT 

g=i 


-i. 

? s ! 


g 


(5) 


and a total-pressure equation 



« P 



(3) 


I r 




R. 
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If minor inert gases are present, Pi^ert " constant, and, therefore, 
Ap^ert ■ 0. The p^_, p 2 , . . ., p g , . . ., Pg. on the left sides of 

equations (10) and (ll) are the species previously selected as inde- 
pendent variables and which appear on the right side of equations (5) . 

Equations (10) , (ll) , (5) , and (3) in the forms indicated may be 
solved simultaneously. Values of Pj_ at temperature T (Pj_^p) may he 

obtained provided that the values of p-j_ at some nearby temperature T 0 
(Pi,T 0 ) are available. 


Knowing *.n the partial pressures (p^rp^) at seme ambient boundary 
where T « T q , select T and guess first approximations for the s in- 
dependent pressures . Using the equilibrium equations, compute first ap- 
pr oximati ons for the remaining Pi^j and then obtain second approxima- 
tions for the independent variables using the Integrated equations (10) 
and (ll) . Iterate to obtain constant values for all Pj^qi which now 


form a new set of boundary conditions. Select a new temperature and 
continue. 


The iteration procedure is normally straightforward and convergent. 
While not always rapid, it can visually be accelerated by successive ap- 
proximations based on the trend of the early steps. For several in- 
tervals in one system that converged in all other intervals, the itera- 
tion did diverge; in those intervals, successive approximations finally 
produced the answer. 

The selection of temperature intervals is not critical. A major 
consideration should be the available tables of Kj_ and Rj_ (there is 

much less work if interpolation can be avoided) . Generally, 400° K in- 
tervals proved satisfactory over most of the range, with 200° K intervals 
near the maximum temperature. It is probably desirable to have the upper 
limit of the h ighest temperature interval within 50° K of the maximum 
temperature. There Is no particular reason for determining the maximum 
temperature. 

At every temperature there are two sets of Pj_ which satisfy the 
equations, one set corresponding to a fuel- lean condition and the other 
to a fuel- rich condition. Accidental interchange of these two condi- 
tions has not occurred; however, a very close approach to the maximum 
temperature might produce such confusion. On either side of the tempera- 
ture maximum It is convenient to arrange equations (10), (ll), and (5) so 
as to use as independent variables those species present in significant 
quantities in that region. 
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To end the calculation for either! the numerical integration or the 
analytical solution , a boundary condition must be selected for the sur- 
face of the fuel particle ; such a criterion may be either the surface 
temperature, the partial pressure. of the fuel, or some equilibrium con- 
dition. The criteria for the work her4 are based on the latter two con- 
siderations. These critieria produced Ipartiele surface temperatures lower 
than those anticipated on the basis of i the Jb oiling point or the sublima- 
tion point. The effects of such differences on the burning rate are rela- 
tively small; however, it is felt that sucii criteria are more realistic 
than efforts to guess the probable surface temperature. 

I T'^ 

Either the numerical integration o{r the use of equations (10) and 
(ll) leads to a plot of partial pressure against temperature for each 
molecular species. Such plots are shown for carbon and boron in figure 
2. The burning rate W f can be obtained frcm-'-a suitable continuity _ 

equation; for example. 


W f " 5? 5 ? CiW± 


( 12 ) 


Substituting the appropriate diffusion Equations (in this case, eq. (2b)) 
gives ’ 


W f 

4jtr^ 



C i 8 i 


dp ± 

dr 


(13) 


Upon integration. 



dr 

4irr^ 



(14) 


for an interval corresponding to a temperature interval AT with limits 
T 0 ' and T, where the partial-pressure limits are p^rp and. Pi,t and - 

the radius limi ts are r»j and rrp. Then ' 
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Summing AT from Tq to T A 



£ °£^Ap ± 


(16) 


corresponds to summing from the ambient boundary rg to the particle 
surface r A . With r c -*■ «>, a value for Wj/4nr^ is obtained.. Summing 

in this way permits the use of an average value of x for each temper- 
ature interval; as a first approximation, * is the thermal conductivity 
of the inert (stagnant) gas from the ambient atmosphere. 


The spatial structure 


w f / w f 
4rtr| 4nr A * 


Then 


of the flame is obtained from the ratio 



4itr A 


fA 

r 


Tg to T 
> jq ' X ? 


Tg to T A 

>_ J * E °iV» t 

AT i 


(17) 


where the denominator is the sum over all AT, and the numerator is suc- 
cessive accumulated sums over AT intervals from the outer boundary to 
the value of the radius r. Flame structures so calculated are shown 
for carbon and boron in figure 3. 

It should be noted that the geometry of the combustion system 
(spheres, c ylind ers, or plates) does not affect the partial-pressure - 
temperature curves involved in these calculations. Geometry is intro- 
duced only in finding the burning rate or the flame structure. Because 
of the mathematical forms of the various geometries, the spherical system 
is most convenient far theoretical treatments. In the spherical case, 
the ambient boundary appears as l/rg, and the term approaches zero as 

the boundary is removed toward infinity; in the other cases, the ambient 
boundary appears explicitly in the final result. In a sense, the steady- 
state solution does not exist except for spheres. 


APPLICATION OF EQUATIONS 
Carbon 

The combustion of a carbon particle in air is treated as a specific 
example of the general equations of the preceding section. The general 
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physical model is shown in figure l(a) , except that no assumptions are 
made about the structure. Carbon, oxygen 1 , and nitrogen are s » 3 atom 
species. C^OjJTq, C^O^Nq, CqOjjNqj CqO^Kq,| and CqOqN 2 are v a 5 molecular 

species, more familiarly, COg, CO, Og, 0,‘ add Ng. The fuel is con- 
sidered as C^OqNq and is not Included as a molecular species because of 

the very low vapor pressure of carbon; tt£s cantssidn merely reduces the 
number of equilibrium equations by one. Then 

w carbon = W C0 2 + W C0 " 1W f * w f 
w oxygen * ^COg + W CO + 2W 0g + w 0 - OW f « 0 

^nitrogen 2 %2 “* ®^f " 0 \ _ 

Numerical integration . - The diffusion equations for v - 1 =* 4 
species, COg, CO, Og, and 0, with the v'tjh species identified as Ng, 

the only inert gas present, take the form. 

W C0 2 - D COg ^COg PCOg^COg d PCOg Deo dp co ^2 ^Og D Q dp 0 \ 

1^2 " ~dr~ * BT~ ~dr~ + W ~~dr~ + RT ’dx~j 

C 2 ^) 

and similarly for CO, Og, and 0. For Ng,' the inert gaB in the ambient 
atmosphere, - 

% " P - £ ;P i C3 c ) 

d i-i : 

The v - s ■ 2 equilibrium equations are written in terms of s * 2 
independent variables (s ■ s minus the mjnnber of inert species; 

3-1*2). Og and COg have been selected dfi independent variables be- 
cause these species predominate at the cold boundary. The chemical 
equations for CO are ; 
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COg ° CO + *2 Og 




CO + 2 °2 “ C0 2 = 0 (zero) 


and for 0 '^molecules" 


°-I°2 


0 - i 0 2 - 0 (zero) J 


(4n) 


The equilibrium equations are 


p 00 2 \ 

Poo - *00 ^72 


p 0 - ) 


(s c ) 


J Z ou 2 

For example, in calculating Fqq at 2000° K using values of KL^ from 
reference 7, equation (6) becomes 


lo S ^co “ log *00 + 2 lo ® ^ 0 % ~ log K C02 

» 21.1878 + j- (6.2695) - 27.1855 « -2.8630 


k - 1.371X10 -3 
CO 


(6q) 


The energy equation is 


w f%„ “ %)%) “ ^CO^O, “ ^2^2 " ^0% " “ “ 4jtr2x < 7 C> 


2 ■ u 2 


dr 


With Wf eliminated from equations (l(j) and. (7 q), the energy equation 
becomes * 


W 0 0 ( B 00 - V + WVi - V + *02% + Vb - *” S ‘ I (V) 
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Use the four ' diffusion equations (2ajj) to eliminate the four s, 
cancel 4nr 2 and dr, introduce 8.j_ e Dj^BT, and rearrange. After 
considerable rearrangement, the energy equation becomes 


[(HCO- 


■ p co " v o) 

H ^ 8 ff22 + f H H \ 8 

' °nn ) °i 


d P0, 


Q 'CO dT tv “C 0 2 “f 0 / “C0 2 die + E 0 2 6 ° 2 dT 


■Vo 


f5a + il. 

dT +i j 


j(Hco -% 0 ) PCO +(HC0 2 “ Hf Q ) Pco 2 +H O e Pp 2 +H ()Poj ^CO ^jr + &0 0 

( 8 C ) 

The continuity equation for oxygen ejan .be similarly treated to yield 


(2 -P, 


CO 


- p o> ( 8 I 


C0 2 dT 0 2 dT 




d PC0 2 ^02, # 

+ B« "JJ + (l + P pri _+P n ) 


L . ^o\ n 

co 2 ‘ *0 2 ' \^C0 “dT“ + 0 Wj 

(19 C ) 


The equilibrium equations are 


v i/ 2 
p ° " ^ ; 

p bo 2 

Pco ■ ^co 


(5 C ) 


"°2 J 

Further, the derivatives of the equilibrium equations are required for 
the numerical integration 


and 


d PC0 . 1 - 

dT " “CO ~TJz 
°2 


£f0 . l/2 ^° " 2 ^2 : 'l dP °2 > \ 

& “ V02 l RT 2 2P0 2 I 




” " B ^ C0 2 . ; 1.. dP 80g ■]_ 

O* : 


ip C 0 2 dT 2p 0 g dT 


apA 

dT I 


? (20c) 
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The energy equation (eq. (18 q)), the continuity equation (eq. (19 q)), 
the equilibrium equations (eqs. (5 q) ) , and the derivatives of the equi- 
librium equations (eqs. (20q)), in the forms indicated, may be integrated 

numerically using the second method of reference 8. This has been done 

liftin g for the cold-boundary conditions (at C) : T = 300° K, 

p 0o ■ 0.2095, Pno 0 “ Poo ** Po " °* 1116 Quantity 5^ (appendix E) was 

assumed to be independent of temperature and composition. Values of 
for equation (6q) and were obtained from reference 7. 

In the numerical integration, the temperature is a very satisfactory 
r unning variable, because and are tabulated at 100° C intervals. 

However, at the imylnturn temperature all dp^/dT - °®. Therefore, the 

running variable must be changed at seme suitably high temperature (in 
this case a shift from T to pg 0 was made, which yielded dp^/dp^o 

and dT/dp C0 ) . 


The boundary condition for ending the calculation of the temperature- 
composition profile was selected as the temperature at which Pqq and 

Pqq 2 satisfy the equilibrium for the reaction 


C + COg 


2C0 I 




p co 

PCOg 


(21c) 


This boundary condition corresponds closely to pqq * 0. The results 

u 

of the numerical integration are values of pgg, PC0g> Po* aQ< l POg a t 

various temperatures (fig. 2(a)); p^ follows directly from equation 
(3q) • 


The burning rate is obtained by substituting the diffusion equations 
into the continuity equation for carbon (the procedure is analogous to 
that yielding equations (12) to (16)): 

"f * 

E PCOo+ECO , ."I 

'C0 d PC0 + ^g^PCOg + — ( 8 C0 d PC0 + ScOg^COg + Sog^POg + 8 0 < 3 ^ > o)J 


C22 0 ) 
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Upon integration, each set of Ap^ corresponding to a particular AT 

corresponds to an Increment in space. Hie sum over all AT intervals 
from Tq to T^ when r c «* » and r^ ! equals the particle radius 

yields the burning rate W f : 


Wf ^ 

4flr A " ^ 


x C B C0AP(X) +8 C0c, A? C0o ) + 


pco 2 +Pco\ 

A % 2 / 


!( B (^PcKD^0q Ap C0o^0o Ap 0o +8 (^P0^ 


( 23 c ) 


p co 2 + p co ; : 

where x and the factor x . are taken as average values over 

-^2 ! 

the temperature interval corresponding to Ap^. The burning rate is re- 
ported in the RESULTS section. 


The complete structure for the carbon particle diffusion flame 
(fig. 3(a)) is traced out in terms of equation (23 c ) as described in 
connection with equation (17) . 


Analytical solution . - To carry out ' the - analytical solution instead 
of the numerical integration, the diffusion equations are taken in the 
form 


W 1 % 4fi 

4 nr 2 " " RT dr 


(2b c ) 


The energy and continuity equations appehr in the form indicated by 
equations (8) and (9) . The appropriate equations can be derived by the 
steps indicated or can be obtained by seating all pi (but not dpi/dT 
or Ap^) equal to zero in equations (18 q), (19q), (22 q), and (23 q ) . 

The energy equation is 

dprn ^CQo ^Oo dp n 

^oCEco - “if" + 8 co 2 ( h co 2 - dr ! + 6 o 2 H °2 "dr" + 8 o B o ar" “ " 1 




tr 


< 8 C> 
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* 


3 

I 


The oxygen continuity equation is 


d 2C0 o 
2B C0 2 dT 


+ 28 o 2 


<^0 2 

dT 


+ %) 


d PCO 

dT 


£po 

dT 


0 


( 9 c) 


The equations may he Integrated at once } provided a smal l enough 
Interval is selected so that the enthalpy terns and 8* may he approxi- 
mated hy average values. On the oxygen-rich side, pg^ and Pco 2 are 

appreciable and are selected as independent variables. The e lim i n ation 
&POg 

of 8 q 2 ■ from the energy equation results in a slight additional 

simplification. Then 


A PCO ? 


AT + Scq^cAPcq + W^Pp 

“ 6 C02 B ifi02 


( 10 C ) 


and 


to 

* 


vhere 


1 /®CO 8g \ 

AP °2 " " APC0 + T Ap 0' + ^CO^COzJ 


. ®Og 

% - =00 - % - -r 




V 


=o ; 


B C0 2 ■ %0g ■ \ - \ 


C Uc) 


Equations (lOg) and (llg) are Iterated together with the equilibrium 
equations (eqs . (5^) ) from the outer boundary toward the region of maxi- 
mum t emp erature, that is, aver the oxygen- rich side of the fl ame front. 
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On the fuel-rich side, pqq and pq^ are the major constituents 

^ 1 dpco 

and are selected as the independent variables. Eliminate S CQ 
from the energy equation. Then ‘ 



at + vs^ + w p o 

- 6 C0 2 H p0 2 


(10c 1 ) 


and 

APc° ■ “ 5^ + 25 02 A P02 + 26 CX)2 Ap C502^ ^O’ ) 

where ; 

*0 " *0 “ *00 + Q 

% m \ m **00 : + ^f Q 
m *X10 z + E £ o ~ ZE CO 

The equilibrium equations (eqs. (5 q) ) must also be rearranged in 

tenns of and p^ 

2 

0 2 - 2C0 2 + 2C0 « 0 (zero) 0|+ CO - C0 2 *» 0 (zero) 

/ p coA 2 soo 2 ^ 5 ° ■* 

p °2 "r° 5oo"/ p °' k *°% ) 


and may be used as determined before but are better evalu- 

ated directly in terms of the Kj_ to avoid a buildup in rounding- off 
errors. This By stem of equations (eqs. (5 q')> (10q')> and (Uq 1 )) inay 
be iterated Just as on the oxygen- rich side _of the maxJjnum temperature. 

The shift from one set of equations ; to the other is best made after 
approaching to within 50° K of the maximum temperature. Using the last 
set of p ± and T found on the oxygen- rich side, guesses are made for 
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S 




■CO 


and p 


COc 


at the same T on the 


procedure is similar, except that with 
small changes in Ap.j_. 


fuel- rich side. The iteration 
AT = 0 it is more sensitive to 


The burning rate is given as 
T c to T A 

4 ^r- ■ 2 w x ( B co^Pco + 8 co ? A Pco ? ) ( 16 c) 

A AS c c 

and the structure of the flame is as in equation (17) . The equations 
can be obtained from equations (22 q) and (23 q) by setting p^ » 0 (but 

not Ap^), except for i ■ Ng» The burning rate is reported in the 
RESULTS section. The structure is plotted in figure 3(a) . 


Boron 

The burning rate of a boron particle is treated in appendix B by 
the methods outlined in the section GENERAL EQUATIONS. Because so many 
equilibrium equations exist, the numerical integration is not carried 
out; rather, the analytical solution using the iterative procedure is 
applied to four separate boundary conditions . The derivation of a 
boundary condition for condensed BgOj is included in appendix B. Also 

included are the details of a rapid three-graph method for obtaining the 
burning rate of boron under a variety of ambient conditions; this 
graphical method is of general use in cases in which some medium temper- 
ature boundary exists. 

Figure 2(b) shows curves of partial pressure against temperature for 
baron burning in air; figure 3(b) gives the flame structure for the same 
case. Hie burning rates from the analytical solution and from the simpli- 
fied equations (appendix C) are presented in the RESULTS section and in 
tables I and. H. 


Hydrocarbons 

Simplified calculations of burning rate have been made previously 
for hydrocarbons in terms of heat conduction and diffusion. References 
1 and 5 give results of calculations showing the effect on burning rate 
of certain changes in the basic assumptions involved. The equations 
used in those calculations can be deduced from the general equations pre- 
sented in earlier sections of this report by introducing certain assump- 
tions about boundary conditions inside the flame. 
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The details of the reduction of the; general equations to the simpli- 
fied equations p re v iously used are contained in appendix C# Ihe use of 
the simplified equations leads to an exp^ssiqn for W f /4jcr A for the 

burning rate, as does use of the extends^. calculation. Ihe only informa- 
tion which wtn be de duc ed about the structure of the flame is the loca- 
tion of the flame front. ! 

'Php burning rate of isooctane (a typical hydrocarbon) in air is in- 
cluded in the RESULTS section. Values for Sther ambient oxygen partial 
pressures and temperatures may be readily obtained from, reference 5 as 
indicated in appendix C. 


Magnesium: 

Work has been published on the experimental burning times (inverse 
rates) of magnesium ribbons in various atmospheres (ref. 6). Ihe calcu- 
lations presented with that work support; the experimental results quite 
closely in the prediction of relative ra^s; however, the calculations 
were based on equations deduced from the : equations used for hydrocarbons 
(refs. 1 5) rather than from the general equations like those pre- 

sented here. If the equations for the magnesium case are deduced from - 
the equations in the GENERAL EOLATIONS section, the resulting equations 
have a somewhat different form than those appearing in reference 6. 

These changes have same effects on the results of the previous calcula- 
tions for magnesium. Hie computed burniijig rates decrease somewhat 
(burning times increase); the predicted Relative rates are somewhat less 
strongly dependent upon the ambient oxygen partial pressure. 

The details of the reduction of the, general equations are contained 
in appendix C. Hie burning rate is obtained as W f /4itr^. The only in- 
formation obtained about the flame structure is the location of the inner 
n.nd outer boundaries of the high- temperature zone . 


RESULTS i • 

When the calculations described in the GENERAL EQUATIONS section 
are applied to carbon as a specific example 7" and to boron (appendix B), 
they yield the burning rate and the flami structure of the diffusion 
flame; these came as simple consequences ; of temperature-ccmpositian 
curves which are the direct result of the numerical integration or the 
iteration of the analytical solution. Temperature-composition curves 
are shown in figure 2(a) for carbon and in figure 2(b) for boron. Ihe 
flame structure of carbon is indicated liji figure 3(a) and of boron in 
figure 3(b). Burning rates are presented in. table I. 


t 
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For carbon, figures 2(a) and 3(a) give the results of both the 
numerical integration and the iterative procedure. They do not differ 
significantly. In the oxygen-rich region below 2100° K, where the dif- 
fusion of 02 and CO 2 is actually equimolar counterdiffusion, the numeri- 
cal integration and the analytical solution (based on diffusion equa- 
tions (2b) of the same form as those for equimolar counterdiffusion) 
give identical plots of partial pressures against temperature. In the 
fuel-rich region, the partial-pressure - temperature curves diverge for 
the two cases j two molecules of CO pass out for each molecule of CO 2 

going in. The burning rates for the two carbon calculations agree well. 
This is an indication that the considerably less difficult analytical 
solution does not introduce significant variations in the results. 

figures 2(h) and 3(b) show the partial-pressure - temperature curve 
and the flame structure for boron burning in air from the analytical so- 
lution (appendix B) . For boron, the analytical solution involving the 
iterative procedure is applied to four separate boundary conditions. 

These results (fig. 4), together with two curves (figs. 5 and 6) repre- 
senting boundary conditions for the condensed oxide, yield boron burning 
rates for a wide variety of ambient temperature and oxygen concentration. 

The following paragraph is an example of the use of figures 4 to 6 - 
finding the burning rate of boron Wj>/ 4*0^ for an ambient temperature of 

400° K in a mixture of 18 percent 0 2 and 82 percent N 2 by volume at 1 at- 
mosphere. From figure 5, the correction term APq 2 « 0.0514; then 

Pq 2 1570 “ °*l so ° “ 0.0514 « 0.1286. Fran, figure 6, for 400° K, 

(Wf/4it) (l/r^yo " l/ r c) a 22x10"® mole per centimeter per second; 
from figure 4, for Pq 2 1570 " O* 1286 stm > (W f /4it) (i/r A - l/r 1570 ) « 
8.33X10”® mole per centimeter per second. Therefore, with l/r^ ■* 0, 
Wf/4jcr£ «* 9.54x10”® mole per centimeter per second. 

The simplified equations deduced in appendix C from the general 
equations are applied to boron. The calculations were carried out for 
two extremely different composition conditions (case A, solely oxygen 
and product, case B, product and nitrogen) ; note remarks on the treat- 
ment of dissociation in appendix C. The values of Wp/ix r^ obtained 

(table i) are in substantial agreement, although the values of Tg 

varied widely. The only significant difference in the two cases arises 
when ambient temperature and oxygen partial pressure are varied 
(table H). 

The burning rates for boron obtained from the extended calculations 
are shown in table U. For a series of increasing ambient temperature, 
the burning rate in air increases roughly 1 to 2 percent per 100° K; 
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reducing the initial oxygen partial pressure from 0.21 to 0.15 atmos- 
phere at 800° K decreases the burning reyte by roughly one- third. Case 
A of the simplified calculation gives numerical answers rather close 
to those of the extended calculation, but there is little variation with 
ambient temperatures and oxygen partial I pres sure. Case B of the simpli- 
fied equations gives numbers about 50 percent greater than those of the 
extended calculation; however, the relative dependence of the burning 
rate upon ambient temperature and oocyge^i partial pressure is similar for 
case B and the extended calculation. When' the results of the extended 
calculations are compared with those of | the - simplified equations, it 
appeal's that for many purposes the simplified calculations should prove 
adequate . 

The numerical results of the calculations described yield burning 
rates which may be expressed in a variety of ways. For the spherical 
geometry considered here, this leads toi values of W f /4tfr A which in- 
clude the steady-state assumption that . l/r c has gone to zero. The 

units of the integrals sire (moles of fuel) /(cm) (sec) . These results axe 
converted to the more physically significant units of fuel consumption, 
g/(cm)(sec), or heat-release rate, cal/(cm) (sec), by introducing the 
molecular weight of the fuel or its moljar heat of combustion. Table I 
gives values of W f /4nr A in a variety of Units for carbon, boron, iso- 
octane, and magnesium. 

The evaporation constant and particle lifetimes for 1-micron par- 
ticles are also included in table I. These values are based on the 
assumption that these particles obey tHe formula for evaporation, which 
Includes burning as a special case (ref-s. 2 and 5), 

d 2 = d| - ifJ't 


where (Iq is the drop diameter at t 4 0, d is the diameter at t, 
and p’ is the evaporation constant. ; The droplet lifetime is then 


tj - 

— M We : 

where B" = 8X10° — t - — square microns per second. 

p 4rtr A 



On a weight and heat-release basib (table i) isooctane and mag- 
nesium have rates which are roughly thtee to five times greater than 
those of boron and carbon. (Molar rates are not a meaningful compari- 
son for substances of different moleeujLar weight.) In teims of particle 
lifetimes, those of boron and carbon ajre an order of magnitude greater 
than those of isooctane and magnesium.' 


* 
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The "burning- rate calculations (which apply only to steady-state 
burning) do not in themselves explain, the extreme difficulty of burning 
elemental boron (e . g. , incomplete combustion even in a bomb calorimeter 
(ref.. 9)). The computed, structure does indicate large concentrations 
of gaseous boric oxide (BgOj) near the boron surface. During ignition 

or under convective conditions , a momentarily cool boron surface might 
become coated, with condensed oxide j such a coating could well preclude 
complete combustion. 


SUMMARY QF RESULTS 

A system with general equations for computing the bur nin g rates of 
anwii 1 particles has been developed and applied to the c ambus tion of 
carbon, boron, magnesium, and isooctane, with the following results: 

1. Burning rates of liquid hydrocarbons and relative rates for mag- 
nesium ribbons have been predicted in satisfactory agreement with experi- 
ment from a model in which heat- and mass - transfer processes are assumed 
rate determining. 

2. Hie equations which determine the steady- state burning of liq uids 
or solids can be set down in a general form. Prom these general equa- 
tions, together with appropriate simplifying assumptions, the explicit 
burning-rate equations which have been used previously for hydrocarbons 
and magnesium may be derived. 

3. If dissociation is considered in the calculations, more realistic 
temperature and composition profiles are obtained for the combustion 
zone. However, such considerations do not change bur n i n g rates greatly 
from those calculated from highly simplified models. 

4. Calculations indicate that heat-release rates are in the order: 
liquid hydrocarbons > magnesium > boron * carbon. 

5. The burning- rate calculations (which apply only to steady-state 
burning) do not in themselves explain the extreme difficulty of burning 
elemental boron. 


Lewis Plight Propulsion laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 14, 1956 
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APPENDIX A 


SYMBOLS 

The following symbols are used in this report (the units indicated 
are those employed here , but the equations are valid for any consistent 
set of units): 

B,B',C definite integrals defined in text, moles/ (cm) (sec) 

C differential function defined by equation (D3) 

c,c 1 constants, developed explicitly in any specific example 

D diffusion coefficient, sq cm/sec 

d diameter, microns i 

j(T) right side of equation (32g). 

G mass -transfer rate per unit £rea, moles/ (sq cm) (sec) 

g general independent moleculaf species 

H enthalpy. Including chemical ■ energy, cal/mole 

AH sensible heat of fuel to boiling point, cal/mole 

H',H" sums of enthalpies, develop e4 explicity in any specific ex-, 

ample, cal/mole 1 

3 * heat of c ambus tion/heat of vaporization 

K equilibrium constant for reference reactions 

k equilibrium constant for specified reaction 

L heat of vaporization, cal/mole 

M molecular weight, g/mole 

n stoichiometric coefficient 

P total pressure, atm ! 7 “ 

p partial pressure, atm 

Q heat of combustion at reference temperature, cal/mole 


r 


s 

sa 
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2 


universal gas constant f 82 ' 05 ( c" om) (atm) /(“O (mole) 

IJL.987 cal/(°C} (mole) 


r 

s 

¥ 


T 

t 

H 

v 

w 

X,Y,Z 


radius from center of particle, cm 
number of atom, species 

s minu s the number of Inert species; number of molecular 
species adopted, as Independent variables 

absolute temperature, °K 

time, sec 

lifetime of particle, sec 
diffusion velocity, cm/sec 
mass -transfer rate, moles/sec 
atom species 

general molecular species 


x 

P' 


mole fraction, p/P 

D^j P/RT, mole/ (cm) (sec) 

evaporation constant, sq microns/Bec 


8 D/xRT, (moles) (°C)/(cal) (atm) 

x thermal conductivity, cal/(sq cm) (sec) (°c/cm) 


v 


number of molecular species 


P 


fuel density, g/ cc 


aUA)*, 

s(2,2)* 


molecular diameter. Angstroms 

Integrals for calculation of diffusion coefficient and 
thermal conductivity, tabulated In ref. 10 


■» 


Subscripts : 

A,B,B',C boundaries of zones (fig. l) 
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ABjBB’ jBCjB'C 

zones between appropriate, boundaries 

(c) 

condensed phase 

eq. 

equilibrium 

f 

fuel 

(g) 

gaseous phase j — 

1,J 

general molecular species 1 

0 

initial conditions 

P 

combustion products 


I* 
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APPENDIX B 


DETAHS OF BORON TREATMENT 

This appendix presents the details of the treatment of boron 
in terms of the general equations . This includes the equations for the 
analytical solution, as well as the derivation of a boundary condition 
for condensed Bg0 3 . Further, there is an explanation of the rapid three- 

£2 graph method for obtaining burning rates of boron under a variety of am- 

§ blent conditions; an example appears in the RESULTS section. 


General Boron System with Gaseous B 2 O 3 

The general physical model is that indicated in figure l(a) , except 
that no assumptions are made about the structure of the flame. Beat and 
mass transfer are considered; convection is assumed to be absent. As 
distinct from carbon, boron is significantly volatile at flame tempera- 
ture. Nevertheless, the handling of the equations is quite analagous to 
the carbon case. 

Boron, oxygen, and nitrogen are s » 3 atom species corresponding 
to X, Y, Z of the general treatment. B, Bg, BO, Bg0 3 , 0, 0 2 , and Ng 

are v = 7 molecular species corresponding to B-^OqNq, B 2 0 qNq, B^O-^Nq, 
BgOjNQ, BqOjNq, Bq0 2 N 0 , and BQOgNg an the notation of the general case. 
BiOqNo ie assumed to be the fuel molecule. 

The continuity equations are 

%oron ** W B + 2 %g + %0 + 2 W BgC >3 “ lw f “ w f | 

w 0xygen * W B0 + ^203 + w 0 + ^2 - 0W f - 0 
^Nitrogen - 2 W N 2 “ CW f - 0 

From the diffusion equation (eq. (2b) ) , the v - 1 ■ 6 diffusion 
equations are of the form 



% Dg dpg 

" “ RT dr 

4itr 


(2bg) 


and similarly for Bg, BO, Bg0 3 , 0, and Og 
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The total-pressure equation for nitrogen is 


P * 2 “ P - 


Pi 


(3b) 


*2: 


For independent -variables, s «* 2 (s m 3 minus 1 inert species, N 2 )$ 
O2 and B2O3 are the noninert species with: greatest partial pressures 
near the cold boundary* The v - s «* 4 Chemical and equilibrium equa- 


tions are 


I B 2 0 3 = B + | 0 2 
3 1 

B + ^ O2 - g 1 ^2^3 = 0 (zero) 

b 2 °3 = b 2 + f °2 

b 2 + \ ~ b 2®3 * ® (zero) 

I b 2 0 3 = BO 4 - i o 2 

50 + \ °2 " \ B2O3 = 0 (zero) 


1 

2 


i 0 2 = 0 


0 - J 0 2 = 0 (zero) 


l/2 

%°3 


Pb * *b “3/4 


^®2°3 


^2 " ^2 1372 


?(%) • 


p o ; 


PBO " fc B0 -jji 


po; 


Po ■ *0*5! 


i/ 2 




r ( 5 b) 




The energy equation tabes the form 


"A. - "A - W ■ “Bcfto - ’Wb„o„ •- w cft> - "oA, - - 4 ”- 2x f 


2 -“2 


2 3 2 3 


2 2 
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As in the general equations, the energy equation and the oxygen continuity 
equation can be transposed to give 


- V ar + ^ ar + 


. v ^*PBoO» dprv 

8 B 2 0 3 ( H B 2 0 3 - m ' + 5 0 H 0 — + 8 0. 




•3-POc 


d PB 0 

dT 


dT 


dT 


2 H °2 dT 


( 8 b) 


and 


^BO 


d PB 0 

dT 


+ 3 




3-PB 2 ° 3 


>3 dT 


dPO 

+ *0 dT - + 28 0. 


*P 0 2 
2 dT 


- 0 


Ob) 


Equations (&q) and (9g) may be integrated using mean values of 5^ and 
HL^, and each equation solved for one of the two independent variables. 

As in the carbon case, it is convenient to eliminate the second in- 

^E^Bg03 

depend^ vartoble fra. tie energy elation, -gj- Is sldsdsated. 

Independent variables are P2g0 3 POg 0111 ^ tie oxyB ® 11 ” 1 *® 11 side, 

Pg^Q^ and Pqq an the fuel- rich side. 


The energy equation (on the oxygen-rich side) is 

AT + BgH^pg + BggS^gApBg + BgO^Sd^PBO + 8 0 H £A p 0 


Ap 0c 




( 10 b) 


The continuity equation is 

APb 2°3 " 35 B 2 o 3 ^ a Bd ip B0 + 5 C^ P 0 + 25 0 2 Ap 0 2 ^ ^1%) 

Equations (IO 3 ) and (llg) together with the equilibrium equations 
(eqs. (5 b)} and the total-pressure equation (eq. (3g)) for nitrogen may be 
iterated, as in the general case and the carbon example, as soon as cold 
boundary conditions are available; however, such a solution is possible 
only when the system is entirely gaseous. 

C o n de nsed B 2 O 3 , a Special Boundary Condition 


When con de nsed B 2 O 3 is present, temperature is relatively low; B, 
B 2 , BO, and 0 are essentially absent; and N 2 , O 2 , B 203 (g), and B 203 (c) 
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are present. The symbols "(g) and (c) with Bg0 3 stand for gaseous and 

condensed (liquid or solid) oxide, respectively. For these conditions, 
the continuity equations assume the f orn^ 


2W B 2 0 3 (g) + ^(c) " w f 
^BgO^g) + 3 % 2 0 3 (c) + m o 2 " 0 






%2 - 9 


J 


(24b) 


The diffusion equations for Bg0 3 (g) and! 0g are 

m ' HP ar 
w o 8 *0g 

w"® a*. 

For nitrogen 

P P °2 ! ^g^Cs) 


> 

J 


(25b) 




The equilibrium partial pressure of BgOg 1 b given by 

lo « p ea,E 2 0 3 -.log K B 2 0 3 (g) - log K B20 3 (c) ( 21 b) 

The energy equation is 

w f*f Q - *k 2 o 3 (g)% 2 o 3 (g) “ w B 2 o 3 (c)% 2 d 3 (c) - " w n 2 % 2 *“ 4jce2x ^ 

(28b) 

Combining equations (24g) , (25 b ) * arui ( 28 B^ yields 

^3^2°sC8) " %2°3( C P# + ^2^°2 “ ” ^JtBT 

- - 1. (29b) 
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Since PB 2 0 3 (g) “ Pec^O^ 


< 3 Peq J B 2 03 B B 2 °3(s) ” ^2°3C c ) s 

dT " dT " PeqjB 2 0 3 ^2 ^ 30 B' 

Substituting for dpggOj/dT in equation (29g) gives the following ex- 
pression for dpQg/dlT: 

dP0 2 

" an? ■ 



Peq J B 2 0 3 ^B 2 Q3 ^^2°5(s) 
BT 2 



2 - 


(»fl) 


A second expression for dpQ /an? is obtained by setting 

u 

*B,0,(c) » 0 in the oxygen continuity equation for the boundary at 

fi w' ' 

which the condensed oxide vanishes, then substituting the diffusion 
equations for ^BgO^g) and W 0g : 

dp °2 5 &B 2 Q 5 d '^ B 2 Q 5 5 8 b 2°5 Pe< b B 2 0 5 ^2°5C c )^ 

dT 2 6 q 2 dT 2 Bq,, Rgj2 


(3%) 


A graphical solution of these two expressions for dpQ /dT in terms of 

H 

Peq,B 203 indicates that the two functions are equal at a temperature of 
1570° K, which corresponds to PB 2 O 3 = 0*004365 atmosphere. This is 
essentially a dew point for B 2 0 3 in the combustion system; it is inde- 


pendent of the ambient temperature and the ambient partial pressure of 

oxygen, but would be affected by any variation of % , 8 0 , or 8 j>q . 

o 2 “ 3 

The condition Wg 2 o 3 (c) "0 as a condition at the dew point is somewhat 


arbitrary; however, no diffusion equation can be written for the con- 
densed phase, and there is no basis for assigning a flow velocity other 
than zero to condensing particles of B 2 0 3 . The condensed particles 


are implicitly assumed to move, in other regions, at the velocity neces- 
sary to maintain the steady state. 
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Ifre partial pressure of oxygen corresponding to T » 1570° K and 
P 3 2 03,1570 is obtained by integrating Equation (31g) graphically to 

yield P02,1570i ' “ 



-^P0 2 : " P0 2 C " P0 2 1570 


(33b) 


where F(T) is the right side of equation (31 B ) • Figure 5 is a plot of 
the integral as a function of (the value of Apo 2 is a correction 
to the ambient partial pressure of oxygen (Pq 2 q) for the departure of 
the ambient temperature from 1570° K). ; 


Burning Rate and Flame Structure 

In the regions where condensed B 2 0 £ exists, the continuity equations 
(24 b) can be combined with the diffusion equation for oxygen (eq. (25 b) ). 
This yields 1 


4 4-0 

Wf - - £ W 0 g - 5 x 8 p 2 4*r — 


(34 B ) 


which when integrated from the ambient boundary (T c , POgC* 8114 r c) bo 

the boundary conditions (1570° K, Po 0 157b> Sb * 1 r 157o) Si ves by equation 
(33b) 2 I - 

^^1570 " ^c)" ” 3 x 8 0 2 ^°2 1570 ‘ P °2 c) 

J r]^70 

1 xF(T) dT 

m> 




(35b) 


'TC 


Figure 6 is a plot of the integral as a function of Tq calculated 
assuming x m The value of the integral in equation (35^) would 

yield the burning rate if r 3^570 could be determined; as it stands the 


integral is a quantity which may be added to the 
to give the burning rate. 


quantity 1 — ' \ 

4 *\ r A r 1570/ 
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With the boundary conditions at 1570° K taken as the cold boundary, 
the ana lytical solution approp ri ate to gaseous conditions may be obtained 
that is, the energy equation (eq. (lOg) ) , the continuity equation (eq. 

(llg) ), the equilibrium equations (eqs. (%}), and the pressure equation 
(eq. (3g)) may be iterated as in the general case and the carbon case. 

Htis solution produces a plot of the partial pressure of each com- 
ponent against temperature, as shown in figure 2(b). The appearance of 
the closed loop in the plot of the partial pressure of boric oxide is 
solely the result of the equilibria and the plot, and has no physical 
significance in itself. 

The burning rate Wf can be obtained using the continuity equation 
for boron by substituting the d if fusion equations and integrating. This 
yields 

v 1570 to T a 

^ ~^ x(& R fip B + ESggdpgg + SbcAPbo + 26 B20^^B203^ 

■ L ' 3 AT 

Cl6 B ) 


WfA 

4jt \ r A 


to which may be added the right side of equation (35g) to obtain 


Vl__jA 
£*Ua re J ' 


w f 

4tcta 


is obtained by letting ~ 


approach zero. 


The structure of the boron flame can be deduced frem the curves of 
the partial pressure, as discussed in connection with equation (17) . 


Graphical Method for Burning Bate 


Because the assumption is made that pg, pg^, and pgQ are equal 
to zero at the 1570° K boundary condition, that pg_Q_ 1570 is can- 

u O’ 

stent, and that the inner boundary is defined by an equilibrium, the 
quantity obtained in equation (l 6 g) is determined by the partial pres- 
sure of oxygen at 1570° K. The analytical solution is carried out for 
a series of such oxygen values. Figure 4 is a graph showing integrated 


values for 


tf f/l 1 \ 
r 1570/ 

\ * 1 s -L.J 


from equation (l 6 g) . The straight-line 


appearance seems coincidental. 
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Figures 4, 5, and 6 permit a rapid, ^valuation of the turning rate 
for wide ranges of ambient partial pressijire of oxygen (see example in 
RESULTS section) . In essence ! 

w 

4rtr^ ** 

„ _ i V w f / i _ A 

" 4=it(r^ “ r 1570 j 4jr^ri5 70 “ t q J 
s g < p o 2 ,i57o) + h ( T c) 

where p 0 1570 is obtained, from Pq q; and. Tq by figure 5. h(T c ) 

6 " u * J 1 . 

represents the quantity plotted. In figurfe 6, and. g(p(v> 1570 ) that In 
figure 4. . 
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APPENDIX C 


REDUCTION OP GMERAL EQUATIONS TO SIMPLIFIED EQUATIONS 

The calculations developed in the GENERAL EQUATIONS section and ap- 
plied to carbon as an example and to boron In appendix B require no as- 
sumptions about any Internal boundary conditions of the flame zone. In- 
deed, the structure of the flame follows as a consequence of the 
assumptions mad** about equilibrium and diffusion. When certain assump- 
tions are made about the flame front. In terms of the models in figure 
1, the general equations can be simplified. The simplification is both 
in the form of and the number of equations involved. The resulting 
simplified equations are those used previously for hydrocarbons (refs. 

1 and 5) and for magnesium (ref. 6). An infinitesimal flame front 
(boundary B, fig. 1(a)), where the concentrations of fuel and oxygen are 
zero, is used for the hydrocarbons and boron. A flame front of finite 
thickness (region BB' , fig. l(b) ) Involving a high- temperature diffusion 
zone Is used for magnesium. In both cases, the fuel surface Is assumed 
to be at its boiling point; Ta ■ Tb.p. . 


may be written 


(ZbH) 


In the form 

Wf [% Q " % + * 0*02 “ V 1 !?] " “ 45trS * H ( 7 H) 

Oxygen and product terms may be written so as to allow for dissoci- 
ation in order to maintain realistic temperatures; the most practical 
assumptions (ref. 5) are atomic homogeneity (a constant ratio among the 
atom species regardless of their chemical arrangement) and chemical 
equilibrium among the molecular species. From the ccxnpositi on- 
temperature information, estimates of the enthalpy- temperature function 
can be obtained. Variations in the enthalpy function in equation (7g) 
produce variations in Tg from equation (40g), but such errors tend to 

cancel out in the graphical integration of equation (41g) , where the 
same enthalpy function appears. 


Hydrocarbons and Boron 
Hie diffusion equation (eq. (2b)) for oxygen 

W 0 2 “DOg d P0g -“Og^f 

4^2 

because Wq 2 => -nQgWf 2 in the over-all process. 

The energy equation (eq. (7)) may be written 
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How if the conditions that pq^ =* p f ;» 0 at the flame-front boundary 

are applied to the energy equation, some simplifications occur. For the 
AB region inside the flame front, there is no net flow W of any compo- 
nent other than fuel. The fuel is assumed to proceed uncracked and un- 
reacted to the flame front; therefore, only the fuel contributes to the 
enthalpy term. The energy equation for the region AB inside the flame 
front becomes ; 

W f [Hf 0 - - W f [AH + L + Hf - - - 4rtr 2 x AB g (37 H ) 

In the region BC outside the flame front, | oxygen diffuses inward, and 
products diffuse out; the energy equation may be written as 

Wf [Hf o + no 2 (Ho 2 - Ho 2 ,c) - n p (H p - %,c) + " n P%,c] rt - 4rtr 2 X]3C || 

( 3 %) 

i — 

Wf[ Q + n 0 2 ^^°2 " “ ^p,c)] B “ 4ltr2>t BG f? (3%) 

Q is the heat of combustion of the fuel particle at its initial 
temperature with oxygen at Tq to give products at To- The terms 
nQgtHog - H Q 2t c) and Up(Bp - Hp^) (deccnjposed into their dissociation 

products where necessary) are a correctioh to the quantity of heat to be 
conducted. In equations (37 jj) and (39g), jthe quantity in brackets is 

the amount of heat that must be conducted j through each temperature re- 
gion in the steady state for each mole of Ifuel consumed. 

In the outer BC zone, the diffusion equation for oxygen (eq. (Zbjj)) 
and the energy equation (eq. (39 h)) yield! 

^2 r P °2,B"° r T B =? dT 

n 0 2 RT)t BC . / dP °2 " J p + n 02 (Ho 2 '- Bo 2 ,c) - n p (H p - H p ,c)J 

•'£o 2 ,c T C 


or 


where 


( 4 %) 

The quantity D Q /x^qT is essentially independent of temperature. With 

the enthalpy function known, the right side of equation (40 jj) can be 
evaluated graphically to determine Tb« 
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The energy equations (eqs. (37 h) and (39 h)) can be integrated in 
their respective zones to yield 


ifr (i i\ _ f * 3 x AB «r f Ts 

"W'v'4 p* ♦ 1 - %,»] %, p * - M v"p<% - M 


- B + C 


( 41 b) 


- B + C 


By using the value of Tg obtained frcm equation (40 h), the two 
integrals on the right side of equation (41g) may be integrated graph- 
ically. Values of x must be estimated because the composition is not 
determined; in the AB zone X£g is appropriately take n as that of the 
fuel; in the BC zone is taken as that of the inert present in the 

ambient atmosphere. Errors made in the estimation of n in the enthalpy 
function and in the selection of x for the evaluation of Tg in equa- 
tion (40 h) tend to ca n cel out whan Wf is evaluated by equation (41 h) • 

As l/rc "*■ 0, equation (41g) yields the burning rate Wf/4*rA. The 
only information easily deduced about structure from this calculation is 
the value of re* which can be determined from Wf anfl the value of 
either of the two integrals. A more compact expression is 


r A 


B + C 


r A B „ 
— - B + C 
r c - - 


B + C 
C 


(42g) 


The form (B + C)/C arises as l/r^ -*■ 0. 

Values for the burning rate of a typical hydrocarbon (isooctane) 
are Included in the EEBULTS section. (They were obtained by use of the 
integrals B and £ frcm figures 8 and 9 and Tg from figure 10 of 

reference 5.) lie values of the integrals in reference 5 are converted 
to the units used in this report, (moles of fuel) /(cm) (sec), when multi- 
plied by 0.1304 (g/lb) (cm/ft)" 1 (g/mole)“ 1 . 


Magnesium 

Hi g h - temperature dissociation data are not available for magnesium 
oxide. Further, there is some evidence that the ood.de dissociates com- 
pletely upon vaporization (this implies that the heat of vaporization, is 
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equal to the heat of reaction at the boiling point). Therefore, for 
magnesium, the flame front of infinitesijmal thickness was replaced (ref. 
6) by a hi gh - temperature zone of diffusion and reaction at the boiling 
point of the cod.de (region BB', fig. 1(a)). 

She . energy equation may be written : 

¥f 0 " W Mg E Mg " “ %0%g0 ■ " 45tr2x |f | 

j 

W f ■ %g + n 0 2 H °2 “ “MgQ^gp] “ - 4ltr2 * ff ) 
Dissociation is considered only in the reaction zone. 

In the AB region inside the reaction zone, there is no net flow of 
any component except fuel. Products are; prtesumed to pass outward (de- 
spite minor oxide deposits on the ribbon; in the experimental case) . Hie 
energy equation in the AB zone becomes > - 

W f [AH + L + Bjjg - § (STfc) 

In the region B'C outside the reaction zbne, oxygen diffuses inward, and 
MgO floats out. The energy equation may be written, analogously to that 
for the hydrocarbon case, as 

Wf [Q + nOgCHOg - Ho 2 , c ) - DMgO(%gO; - %gO,c)]- - 4xr 2 *B>c |jf 

(39jjg) 

Because of the condensed product, the diffusion equation for oxygen 
reduces to the equation for the diffusion of oxygen through a stagnant 
film of the ambient inert gas (form of eq. (D6) , appendix D) : 

w o 2 -Bo 2 f , «f0 2 

tar 2 ■ unr - Pb 2 >. ^ 


(«Mg> 
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The diffusion equations for oxygen and the energy equation for the 
outer B'C zone, combined and integrated, yield 


^2 

n 0 2 ETx B 


£ 0 2 ,B-=? Pdp02 , 

'C I P-Po 2 J 

‘ / P0 2 ,C ^ 


^?B ,S3 ^b.p. 


dT 


|Q+n0 2 (Ho 2 -Ho 2 ,c) 

(^^Mg) 


In this case, the graphical integration yields the partial pressure of 
oxygen at the B* boundary rather than the temperature of the flame front, 
as in the hydrocarbon case (eq. (4 Qh)). Ebg/^CB'C again assumed to 

be independent of temperature. 

In the high- temperature diffusion zone (EB 1 ) combining the con- 
tinuity equations 


2W0g + WjjgO =* 0 
%g + W MgO " W f 


(^-Mg) 


with the energy equation (7^) yields 

V%0 ■ H Ms ) + 2W ° 2 (%gO - • \ H ° 2 ) ” ‘ 4,tA § 

In the BB' region of high-temperature diffusion, the temperature gradient 
is zero; 


* - - ^ 




and similarly 


( 46 Mg) 


W f - - 2(J^+ 1)W Q , 




Here 03200 Indicates the heat of combustion of magnesium at T 0 to' 
form condensed MgO at 3200° K. The factor “ 2.0 does not appear in 
reference 6. It may be regarded as the ratio of the amounts of MgO con- 
densed at the outer and inner boundaries of the BB 1 zone. 
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In the high- temperature region, the 1 diffusion equations for oxygen 
and magnesium are in the form of equations (2a) . Both diffusion equa- 
tions and both forms of equation (46^) combine to yield 

W f D °2 F (2 2) <3p °2 Ul v 

4jtr 2 " RT P + (2>;- l)Po 2 ^ 

Equation (4?Mg) is integrated' over the boundary conditions po 2 ■ P0 2 ,B' 
at r * rgi and p 0 _ » 0 at r ■ r-g, 'In Reference 6, the factor 

o 1 

(2Jf" + - l) * 2.0 does not appear in the integrated form of the 

diffusion equation; rather the factor ib l/n^ ■ 2. The introduction 

of from the general equations, instead of the use of equation (43jjg) 

for the diffusion of oxygen in the high- temperature zone, does result in 
some modification in the dependence of the burning rate upon ambient 
partial pressure of oxygen; the original formulation, which now must be 
regarded as semi empirical, fits the experimental results somewhat more 
satisfactorily. 

Finally, the burning rate is obtained by adding the integrated forms 
of equations (37^) , (39^) , and (47^) . As l/r c 0, 



Each term of equation (48^) is positive with the logarithmic term domi- 
nating, particularly in oxygen-rich atmospheres. Equation (44jjg) shows 
a strong dependence of Po 2# b' u P on Pp 2 ,C> £un ^ ) i en ^ partial pres- 
sure of oxygen; this shows up in turn ip. equation (48^) as a strong de- 
pendence of Wf upon Pog,B' * is i consistent with the experimental 

evidence for ambient atmospheres corverijog a range of 15 to 100 percent 
02 (ref. 6). 
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APPENDIX D 


4l 

o 


DERIVATION OP APPROXIMATE DIFFUSION EQUATIONS FOR 
MULTICOMPONENT MIXTURES 

General equations for diffusion in multicomponent mixtures are pre- 
sented in reference 10 (pp. 516-520) . If external forces, pressure 
gradients, and thermal diffusion are neglected, equations (8.1-3) of 
reference 10 became: 


E SSrCVy 


J-l 





1 3 2j m • m j V ■ 1 


(Dl) 


where n^ and nj are concentrations of species i and J in mole- 
cules per cubic centimeter ;_n is the total number of molecules per 
cubic centimeter; Vj ' and are diffusion velocities, in centimeters 

per second; D^ is the binary diffusion coefficient; and r is dis- 
tance in centimeters. In a v component mixture, only v - l of these 
equations are Independent. 

The diffusion velocity V-^ or Vj is related to the flux of mole- 
cules G-^ or Gj (moles /(sq cm) (sec)) through a unit area normal to the 
direction of diffusion. Assume the gases obey the ideal gas law; then 



RT 
xj P 


where xj =» nj/n is the mole fraction of component j. A similar equa- 
tion relates V^ and Gj,* 


Equations (Dl) became 

*i E i£r - °i E 


J“1 


J-l 




£fi 

dr 


1 , 2 , 


v - 1 


(D2) 


where = D^P/RT* If, as an approximation, it is assumed that the 
binary diffusion coefficients of component i into each of the other 
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components of tbs mixture are equal, then 
* • • ** P j_ ^ ® Pjl_ , and 




*1 ^ G i “ G i 1 -It 


or 


V-l 

% Z) g j 


dxi 

- «i » P ± - *A 


■ p iJ = 


; i * 1,2, . . ., v - 1 (D3) 


where G v is the flow of the component :for which no equation has been 
written. 


It is convenient to solve these equations (by means of determinants ) 
for GjL as a function of all the xj a£d Cj : 



(D4) 

Obese equations are strictly valid only lin mixtures where binary diffu- 
sion coefficients between all the pairs 'of species are equal, that is, 
in binary mixtures. 


i 
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When there is an inert gas present in the system, it is stagnant 
(its flew is zero) . In problems concerning the burning of liquid or 
solid fuels in air, nitrogen is such a diluent, provided formation of 
nitrogen atoms and NO molecules is neglected. In this event it is con- 
venient to select G v = Gw o 0. (This is not convenient if the amount 

u 

of inert gas is small, because as p v -+ 0,pj_/p v -*■••) Equation (D4) 
then becomes 


-°i 


Di dpi Pi ^ Dj d Pj 

KT dr + Pw 2—d RT dr 
n 2 


(D5) 


In a binary mixture, pw 

u 


P - Pi, and equation (D5) becomes 


-G ± 



api 
RT dr 


(D6) 


which is the familiar equation for diffusion of one gas through a stag- 
nant film, of a second gas (ref. ll) . 


Finally, when the concentrations of diffusing gases are small rela- 
tive to the amount of inert gas, the second term may be neglected 
(pi/p v -*■ 0) . In this event the simplest form of diffusion equation 
occurs : 


„ &L ^i 
“ Cx i " HE dr 


(D7) 


which is the same fona as the equation for equimolar counter diffusion 
in a binary mixture. 
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APPENDIX E 


ESTIMATION OF \ FCR CARBON AND BORON COMBUSTION 


The quantities 5^ are defined by 


Dj. 

81 E xRT 


(El) 


where D* is the diffusion coefficient for component i into the gas 
mixture, and x is the thermal conductivity. For the combustion of 
carbon and boron in air the properties of the gas mixture exclusive of 
component i are approximated by those of nitrogen. From the expres- 
sions for diffusion coefficient and thermal, conductivity presented in 
reference 10, 



where A 


i-Nc 


ffl (2,2)*/ a (l,. 
1-N 2 / i-N 


l)* 

2 


The ratio 2 


( 2 , 2 ) 


N, 


/“i-N; 


2) H 


is further approximated as equal to 


unity. The quantity A£_ N2 is close to 1.1 for a Lennar d-Jones (6-12) 


interaction potential between molecules (ref-. 10 ) . The heat capacity of 
nitrogen has been taken at 2000° K (ref. 7) (a reasonable mean tempera- 
ture) . With these approximations 



In order to estimate &j_ values, the molecular diameters for the 
atoms and molecules involved must be obtained. 
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For carton combustion, a values for Ng, Og, CO, and COg were taken 

from reference 10. For atomic oxygen, the diameter was estimated assum- 
ing the same mean density of matter in the atom as in the molecule, that 
is, 


atom 


/I 

f 2 diatomic molecule 


0.794 a 


molecule 


(E4) 


(Thi s method probably underestimates somewhat the diameters of atoms, 
since the electron density of an atom is usually lower than that of the 
correspon ding diatomic molecule.) Diameters chosen and the correspond- 
ing 8 values at atmospheric pressure are shown in the following table: 


Molecule 

Diameter, 

Diameter 

8 at 1 


O 

obtained 

atmosphere 


A 

from - 

(eq. (E3)) 

% 

3.749 

1 

0.1192 

°2 

3.541 

fflef. 10 

.1220 

CO 

3.706 


.1205 

C0 2 

3.897 

) 

.1037 

0 

2.810 

Eq. (E4) 

■ 1827 


For boron combustion, diameters of Eg, Og, Bg, BO, 0, B, and Bg0 3 

must be obtained. (Literature values ex±Bt only for nitrogen and 
oxygen.) It was found, however, that a plot of the quotient of molecu- 
lar diameter (ref. 10 ) and intemuclear dista n ce (ref. 12 ) against in- 
temuclear distance for a number of diatomic molecules gave a smooth 
correlation. This is shown in figure 7. From this correlation and the 
known internuclear distances of Ng, Og, Bg, and, BO (ref. 12) molecular 

diameters for the diatomic molecules were obtained. Atomic diameters 
for B and 0 were then calculated from equation (E4) . The diameter of 
B 2 O 3 was estimated from the density of crystalline B 2 O 3 , with the mole- 
cules assumed to occupy 74 percent of the crystal volume (closest pack- 
ing of spheres) . Results are shown in the following table: 


Molecule 

Diameter, 

1 

Diameter obtained from - 

8 at 1 atmos- 
phere (eq. (E3) ) 

&2 

3.70 

) 

0.1192 

02 

3.52 

(Fig. 7 and internuclear 

.1212 

B 2 

3.82 

| distances of ref. 12 

.1238 

BO 

3.50 

J 

.1273 

0 

2.80 

}ecl. (E4) 

.1812 

B 

3.04 

.1925 

B 2°3 

4.50 

Density of B 2 O 3 (crystal) = 
1.805, and assumption that 
74 percent of volume is 
occupied 

.0812 
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TABLE I. - COMPUTED BURNING RATES OF VARIOUS SUBSTANCES 


[Temperature, 300° K; air at 1 atmosphere.] 



Burning rate 

, W f 

Evaporation 
constant, P', 
sq. microns 

Lifetime of 

1-micron 

particle. 


moles 

gm 

cal 


(cm) (sec) 

(cm) (sec) 

(cm) (sec) 

sec 

*19 

sec 

Isooctane a 

0.86X10-5 

99X10-5 

10.5 

1.11X10 6 

0.90X10“® 

Magnesium 8- 

3.97 

97 

5.7 

AAA 
• HU 

2.25 

Boron 
a Case A 

1.31 

14.2 

1.98 

.0492 

20.3 

a Case B 

1.79 

19.4 

2.70 

.0673 

14.8 

Analytical 

1.27 

13.7 

1.92 

.0477 

21.0 

Carbon 

Numerical 

1.05 

12.7 

1.02 

.0505 

19.8 

Analytical 

1.10 

13.2 

1.07 

.053 

18.9 


Simplified calculation. 


table: n. - effects qf ambient temperature and 

AMBIENT PARTIAL PRESSURE OF OXYGEN 
ON BURNING RATE OF BORON 


[Pressure, 1 atmosphere. ] 


Temper- 

ature, 

°K 

Oxygen con- 
centration, 
percent 

Burning rate, 
moles/ (cm) (sec) 

Extended 

calculation 

Simplified calcula- 
tion (eq.. (41g)) 

Case A 

Case B 

300 

20.95 

12.7X10-6 

13.1X10-6 

17.9X10-6 

800 

20.95 

13.2 

13.1 

19.1 

1300 

20.95 

14.5 

13.0 

20.4 

800 

15.00 

8.2 

11.2 

12.1 





Ambient atmosphere 



(») 


figure 1. - Crass- sectional models 


Ambient atmosphere 



burning fuel particles. 
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co 

& 



Figure 2. - Pres sure- temperature ourvea for burning in air. 
Arrases Indicate outward direction. 


Partial pressure, p. 
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Temperature; °K 
(b) Bor oil. 


Figure 2. - Concluded. PresBure-temperaturie curves for burning In air. Arrows 
indicate outward direction. ' 
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Ambient teiverature^ °K 


Figure S. - Correction to ambient partial pressure of oxygen for 
temperature from 1570° K (eq. (54^) , appendix B) . 
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Molecular diameter 

From viscosity 
From second virlal 
coefficients 


trogen 


'bon TnoTuojd.de 


ltrlc oxla 


parogen] 

tilorldeJ 


3gen Iodide 


-Chldrine-I 


1.0 1.2 


1.4 1.6 1.8 Q 2.0 

Ibternuclear distance, A 


2.4 2.6 


Figure 7. - Correlation between molecular diameter (ref. 10) and lntemnclear distance (ref. 12) for 
diatomic molecules . 
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